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The unified model of active galactic nuclei (AGN) predicts silicate emission 
features at 10 and 18 /zm in type 1 AGN, and such features have now been ob- 
served in objects ranging from distant QSOs to nearby LINERs. More surprising, 
however, is the detection of silicate emission in a few type 2 AGN. By combining 
Gemini and Spitzer mid-infrared imaging and spectroscopy of NGC 2110, the 
closest known Seyfert 2 galaxy with silicate emission features, we can constrain 
the location of the silicate emitting region to within 32 pc of the nucleus. This is 
the strongest constraint yet on the size of the silicate emitting region in a Seyfert 
galaxy of any type. While this result is consistent with a narrow line region 
origin for the emission, comparison with clumpy torus models demonstrates that 
emission from an edge-on torus can also explain the silicate emission features and 
2-20 jum spectral energy distribution of this object. In many of the best-fitting 
models the torus has only a small number of clouds along the line of sight, and 
does not extend far above the equatorial plane. Extended silicate-emitting re- 
gions may well be present in AGN, but this work establishes that emission from 
the torus itself is also a viable option for the origin of silicate emission features 
in active galaxies of both type 1 and type 2. 

Subject headings: galaxies: Seyfert — galaxies: individual (NGC 2110) — in- 
frared: galaxies 



1. Introduction 



In the unified model of active galactic nuclei (AGN), the differences between objects of 
types 1 and 2 are explained by a torus of dust and gas that obscures the broad line gas from 
some viewing angles while leaving it exposed from others. Variou s lines of evidence, notabl y 
the detection of polarized broad emission lines in type 2 objects (lAntonucci &: Miller! Il985l ). 
support the hypothesis of anisotropic obscuration, but the chemistry, structure and origin of 
the torus, as well as the true extent of its role in unifying the various types of AGN, remain 
the subject of intense study. 

The dust in the torus is expected to absorb short wavelength radiation from the nucleus 
and re-emit it in the mid-infrared (MIR), and a strong constraint on models aiming to 
explain and predict emission from the torus is the Si-0 bond stretch near 10 /iin. In the 
simplest realisations of the unified model, an edge-on view through cool dust in the torus 
causes a prominent absorption feature in type 2 AGN. In type 1 AGN, where hot dust 
at t he inner surface of the torus is revealed, models often predict strong silicate emission 
(e.g. lEfstathiou fc Rowan-Robinsonlll995l : iGranato Danesdll994l ). depending on the torus 
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geometry (IPier fc Kroliklll992l ). 



In practice, however, the situation ha s proved less straig htforward than this. For in- 
stance, in a sample of 16 Seyfert 1 galaxies, iRoche et al.l (ll99ll ) detected silicate emission in 
only one. More recently, Spitzer observations of Q SOs have establ i shed not only that the sili 



cate emission feature is common in those objects (IHao et al.ll2005l . 120071 ; ISiebenmorgen et al. 
20051 ) but also - surprisingly - that silicate em ission is present in the average type 2 QSO 



spectrum ( jSturm et al.l I2006T ) . IShi et al.l (120061 ) find that in a Spitzer/IRS sample of ~100 
AGN, X-ray column density roughly correlates with silicate feature strength, in agreement 
with unified schemes. However, there are some obvious anomalies: type 1 AGN with strong 
silicate absorption and type 2s with silicate emission. These observations are not readily 
explained by simple torus models. 



Sturm et al.l (120051 ) note that the dust temperature implied by the silicate bands in the 
QSO spectra is ~200 K, suggestive of dust in the narrow- line region (NLR) rather than 
in the hot inner wall of the torus. NLR dust with a temperature of ~200-300 K has been 
dete cted in mid-IR images of several AGN, extended over ~l-2" from the unresolved nucleus 



fe.g. lBock et alJl200d iRadomski et al.ll2003l : iGoriian et al.ll2004 IPackham et al.ll2005l ). Such 
dust, heated by the central engine but not directly associated with the torus itself, could ex- 
plain the unexpected silicate emission in type 2 objects. Alternatively, alterations to the basic 
torus models may enable them to match t he observations. Several recent torus models have 
incorporated a clumpy dust distribution dNenkova et al.l |2002| ; iDullemond &: van Bemmel 
20051 ; iHonig et al.l 120061 ; ISchartmann et al.ll2008l ). which is expected from c onsiderations of 
grain survival in the circumnuclear environment (IKrolik fc Begelmanl Il988l ) and whose ex- 



isten ce is suggested by interferometric observations of the Circinus galaxy ([Tristram et al. 



2007). In a clumpy torus, cool dust may be present close to the central engine. iNenkova et al. 



( ]2008bl ) find that in certain circumstances a clumpy torus can produce silicate emission fea- 
tures even for lines of sight close to the equatorial plane, especially if the optical depth per 
cloud Ty > 100 or if the number of clouds along the line of sight is as low as Nq ~ 2. The 
silicate emission in type 2 AGN could therefore come from either the NLR or the torus. 

While Spitzer's excellent sensitivity is well suited to discovering silicate emission in 
AGN, its spatial resolution cannot pinpoint the location of the emitting dust. Ground-based 
mid-IR spectra of NGC 1068 show that spatial resolution is critical: the depth and profile of 



the silicate feature y ary considerably on sub-arcsecond (<70 pc) scales (IMason et al.l 12006 



Rhee fc Larkmll2006l ). Beyond that work, the spectral properties of the s mall-scale extended 



MIR emission have bee n explored in only a handful of Seyfert galaxies (IRoche et al.l 12006 



20071 : lYoung et al.ll2007l ). A detailed dissection of the nuclear regions is necessary to establish 
how — if at all — these results fit within AGN unified schemes. 
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At a distance of 32 Mp c (ff n = 72 km sec" 1 Mpc" 1 ; 1 // =150 pc ), NGC 2110 is the 
neare st known type 2 AGN faradt et al.lll978l : iMcClintock et al"1ll979h with silicate emis- 
sion (IShi et al.ll2006l ). and is unique in being both bright enough for ground-based mid-IR 
spectroscopy, and close enough for spatial resolution of a few x 10 parsecs to be attainable. 
To investigate the nature and location of the region(s) giving rise to silicate emission in 
active galaxies, we present ground-based and Spitzer mid-IR imaging and spectroscopy of 
the nucleus of NGC 2110. 



2. Observations 



N band spectroscopy and an imaging observatio n in the 11.2 /im N' filter (AA = 
2.4 /im, 50% cut-on/off) were acquired using Michelle (IGlasse et al.lll997l ). the mid-IR im- 
ager/spectrograph on the Gemini North telescope, on 20070318 and 2OO7O32C0. Conditions 
were clear and dry (PWV < 1 mm) on both nights and a standard chop-nod observing scheme 
was employed for both imaging and spectroscopy. The 2-pixel (0.36") slit was used, giving 
spectral resolution A/AA ~ 200 and spectral coverage of the whole N band atmospheric 
window. T he slit was oriented at 160°, appro ximately along the inner ionization cones of 
the galaxy (jPoggd Il989l : iMulchaey et al.lll994l ). Because of an instrumental problem at the 
time of the observation, the galaxy was not optimally centered in the slit. Nonetheless, an 
adequate signal-to-noise ratio was achieved and the pointlike nature of the nucleus (below) 
means that the conclusions of this work are not affected. On-source exposure times were 
1200 sec for the spectroscopy and 150 sec for the imaging. 

The galaxy appears pointlike in the short imaging observation, with FWHM^0".42 
(comparable to the FWHM of the photometric standard star) and radial profile well fit 
with a Moffat profile. There is no evidence of significant low-level extended emission in 
the smoothed data, implying that the Michelle spectrum represents the unresolved nuclear 
source. The total flux density of the nuclear point source was measured to be 286 mJy; 
variable atmospheric transmission in the mid-IR means that this value is likely accurate to 
~10%. The Michelle photometry agrees with the IRS spectrum (Fig. [TJ Shi et al. 2006; 
Gorjian et al. 2009, in prep) to within the uncertainties^. 

Initial IRAF processing of the Michelle spectra involved combination of the chop- and 
nod-subtracted data, and use of a median subtraction algorithm to remove channel-channel 



Program ID: GN-2007A-Q-49 

2 IRS spectroscopic flux calibration is accurate to ~ 5-10%; IRS Data Handbook, version 3.1 
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Rest wavelength (jum) 

Fig. 1.— Spitzer/IRS and Gemini/Michelle spectra of NGC 2110 (3.7-10.7" and 0.36" slits, 
respectively). The hatched area shows the telluric O3 band, which is not well removed by 
division by the standard star. The Michelle spectrum was calibrated using the N' imaging 
but has been scaled slightly for ease of comparison with the IRS data. 

offsets and low-level vertical striping. The spectrum was optimally extracted using the 
Starlink Figaro package and wavelength calibrated using sky emission lines in the raw frames. 
The galaxy spectrum was then divided by the spectrum of a G8III star (HR1784) to cancel 
telluric absorption lines, multiplied by a 4960 K blackbody curve, and flux calibrated using 
the imaging data. 



3. Results and discussion 

The Michelle and IRS data on NGC 2110 are presented in Fig. [TJ with continuum- 
subtracted spectra in Fig. [2J The continuum shown in Fig. [T] is a spline fit to the IRS 
spectrum at 5 - 9 jttm and 24 - 30 /im. Although fitting only to the extremes of the spectrum 
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may exaggerate the strength of the silicate emission (jSirocky et al.ll2008l ). we use this fit for 
ease of comparison with previous work. The same continuum was also subtracted from the 
Michelle spectrum after a slight scaling to match the flux density in that spectrum. 

Outside the 11.3 fim PAH band, the 8-13 /xm spectral region changes very little between 
the Michelle and IRS spectra; the flux densities of the continuum and silicate emission in 
the Michelle and IRS spectra agree to within the uncertainties O- The similarity of the 
strength and profile of the 10 /im silicate feature in both spectra implies that the silicate 
emission in NGC 2110 arises within the point source detected in the ground-based imaging, 
whose FWHM < 63 pc. At that distance the emissio n could arise in t he torus itself, whose 
outer radius is thought to be no more than a few pc (jJaffe et al.l 12004 ; iPackham et al.l 12005 



Tristram et al.ll2007l). or in the inner par t of the narrow line region, which extends to r ~ 230 
pc in the optical (IMulchaey et al.lll994l ). The compactness of the emission argues against 



an or igin in mass-losing evolved stars, as detected in Virgo cluster galaxies (IBressan et al. 



200J). 



Fig. [2] also sho ws the silicate features in the LINER/type 1 Seyfert galaxy, NGC 3998 
( iSturm et al.ll2005l ). The features in NGC 3998 are weaker than those in NGC 2110 but the 
profiles and peak wavelengths of their 10 /im bands are quite similar. However, the 18 /im 
feature in NGC 2110 is weaker relative to the 10 /im band than in NGC 3998, and NGC 2110 
exhibits a more pronounced red wi ng on the 18 nm fe ature. In this respect, NGC 2110 bears 
more resemblance to the QSOs in lHao et al.l (120051 ) . which also have strong red wings to 
their 18 /im bands. The peak height of the 18 /im feature is sensitive to the choice of 
continuum, but other reasonable continuum fits give similar 10 /im profiles and also show 
the red wing on the 18 /im band. The properties (profiles, peak wavelengths and relative 
strengths) of the silicate bands in NGC 2110 lie w ithin the range exhibited by other AGN 
with silicate emission features. ISturm et al.l (120051 ) note that the 10 /zm band in NGC 3998 
is broadened and shifted to longer wavelengths than in the Galactic ISM, and discuss grain 
size and composition effects that could explain the observed spectrum. The similarity of the 
10 /im band in NGC 2110 to that in NGC 3998 suggests that such grain processing may also 
be at work in this Seyfert 2 nucleus. 



Schweitzer et al.l (120081 ) have modeled silicate emission arising from dust in the NLR of 



PG QSOs. The distance to the dust giving rise to the feature depends on various model 
parameters (e.g. NLR density, extinction), but on average -Rdust ~ 80(Lboi46) 1,/2 pc, where 



Lboi — 7L(5100) and is in units of 10 46 erg s 1 
from the nucleus would require L hn i 



> 1.4 x 10 4 5 erg sec" 1 



For the silicate-emitting dust to be > 32 pc 
No 5 100 A continuum luminosity 
is available for NGC 2110, but iMoran et al.l (120071 ) find L bol ^ 2 x 10 44 erg sec" 1 based on 
the extinction-corrected [OIII] A5007 line. This implies -Rdust < H pc, consistent with the 
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NGC 2110 

NGC 3998 * 1.6 




10 



20 



Rest wavelength (/xm) 

Fig. 2. — Continuum-subtracted Michelle and IRS spectra of NGC 2110 (in set), and 
continuum-subtracted IRS spectra of NGC 2110 and NGC 3998 (ISturm et al.ll2005l . emission 
lines removed). The NGC 3998 spectrum has been multiplied by 1.6 and the data in the 
main figure displayed on a logarithmic wavelength scale to emphasize the feature profiles. 



constraint on the distance to the silicate emitting dust imposed by the comparison of the 
Michelle and IRS data. This distance is clearly much smaller than the optical size of the 
NLR in NGC 2110, an observation also noted for their PG QSO sample by Schweitzer et al., 
who suggest that the silicate-emitting dust may be associated with highly ionised gas in the 
innermost regions of the NLR. 



The calculations of ISchweitzer et al.l (120081 ) demonstrate that the silicate emission in 
AGN may arise in the NLR, but they do not rule out that some or all of it could be produced 
by dust in the torus (and there need not in fact be a physical difference bet ween the outer 
edge of a clumpy torus and the beginning of the NLR). For clumpy tori, iNenkova et al. 
(]2008bl ) find that if the optical depth per cloud ry > 100 or if the number of clouds along 
the line of sight is iVo ~ 2, a silicate feature can be produced even for lines of sight close to 
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the equatorial plane. 



We fit the clumpy torus models of lNenkova et al.l (j2008al ) to the MIR spectra of NGC 2110. 
In these radiative transfer computations, the individual clouds have a fixed optical depth. 
The additional parameters detail the distribution of the clouds. They are radially distributed 
according to a power law, oc r~ q , from the dust sublimation radius, Rd to an outer radius R Q , 
which we parameterize with Y = R /Rd. The torus has width parameter a with a Gaussian 
edge. The average number of clouds along an equatorial ray is N . The number along a 
given line of sight declines with altitude, with Ni os (f3) = N exp(—j3 2 /a 2 ) along angle (3 from 
the equator, and the emergent reprocessed emission is calculated for all viewing angles, i. 

We fit only the MIR continuum, using the relatively line-free Michelle spectrum from 8 
to 12.5 /im (except the ozone region) and the I RS spectrum in regions ind icated in Fig. [3J 
Because of uncertainties in starlight su btraction (lAlonso-Herrero et al. Il996h and the possible 
variability of this object in the NIR Lawrence et alJll985h . we do not fit the photometric 
data points. Neither do we fit to the LL spectrum of NGC 2110; although most of the 
11.2 fim flux arises in a source with FWHM<0.42", contributions from extended emission 
may become significant at longer wavelengths. 

We allow all the model parameters to be free. One disadvantage of the models is that 
they are highly degenerate. Concentrating on the MIR, where only the clouds close to 
the AGN (within ~ 15i?d) are relevant, a number of different models produce similar fits 
and spectral shapes. The total cloud distribution determines the MIR emission, so various 
combinations of total number of clouds, outer extent, and radial profile produce similar 
spectra. The best-fitting models all have near-equatorial views through the torus, which 
has a small width parameter. The inclination of the torus is driven by the NIR/MIR ratio 
in the data; more face-on models cannot fit the short end of the IRS spectrum without 
predicting less 1 4 jum emission than is observed. Larger tori may produce a slightly lower 



NIR/MIR rat io flNenkova et al 



2008 bJ) but large to ri are inconsistent with high resolution 



observations (jJaffe et al.ll2004 iTristram et al.l 120071 . although see Kishimoto et al. 2009 for 
contrary evidence). 

We plot the formally best-fitting model ("model 1") in Figure [31 which has ry = 20, 
Y = 30, q = 0, N = 5, a = 15°, and % = 80°. A similar model ("model 2") changes N a to 
10 and i to 70°. The tori have smaller opening angles than required by the unified model to 
explain the relative numbers of types 1/2 AGN. However, slender tori are suggested by th e 
disk-like broad line region (BLR) which is detected in polarized light (IMoran et al.ll2007l ). 
The apparent scarcity of type 2 AGN with silicate emission may also indicate that these 
objects are outliers in terms of torus properties. A somewhat different solution ("model 3") 
has fewer clouds along a given line of sight, and the torus covers a larger angular extent 



- 9- 



(r y = 40, Y = 30, q = 0, N = 2, a = 45°, and i = 90°). All these models obscure the 
central engine, with Nh ~ 1-2 x 10 23 cm -2 along the line of sight, assuming a standard 
gas-to-dust ratio. This Nm is within a factor of a fe w of that estimated by most authors 



f lMalaguti et al.lll999l : iRisalitil 12002 ; Evans et al.l 120071 ) 



The vertical scaling is a free parameter and sets the bolometric luminosity of the model 
source. We find L hoi = 4, 3, and 0.7 xl0 44 ergs _1 in models 1, 2, and 3, respectively. These 
yalues are roughly consistent with the L^ \ m 2 x 10 44 erg sec -1 estimated from [OIII] A5007 



( jMoran et al.l 120071 ). In all cases, the silicate-emitting region is small. For a standard AGN 
heating spectrum and a dust sublimation temperature of 1500 K, Rd = 0.4(Lfe o //10 45 ) 1 / 2 pc. 
Thus, the outer extent of the torus we model, at 30 Rd, is < 8 pc. 

One problem with all these models is that they fail to capture the steep rise of the silicate 
feature longward of 10 /im. The wavelength of the silicate peak is a fun ction of the dust 



compo sition. All the models we present use the cold astronomical silicate of lOssenkopf et al. 



(119921 ). which has a maximum optical depth at 10.0 /im. We do not experiment with alternate 
dust compositions, but we do demonstrate in Figure H] that even an obscured AGN may 
exhibit strong silicate emission. "Model 4" is similar to models 1 and 2, with a slender torus 
(<r = 15°). The primary difference is that it contains fewer clouds along radial rays, but 
with nearly equatorial viewing (i = 80°) and ry = 40, the average line-of-sight optical depth 
is still 1 x 10 23 cm -2 . Although this model is formally a poor fit, it has the advantage that 
the shape of the 10 /im silicate feature does not show the "double-peaked," self-absorbed 
structure seen in models 1-3. The dust composition determines the shape and strength of 
the silicate feature, while variations of the grain size affect the relative extinction of the 
continuum. A distribution favoring large grains does not alter these models significantly. 
Favoring small grains corresponds to smaller values of Ty for the same silicate strength and 
profile. The NIR extinction is then reduced relative to the MIR, so the preference for edge-on 
views remains robust. 

To summarize, these observations permit us for the first time to set tight limits (r<32 
pc) on the size of the silicate-emitting region in a Seyfert galaxy. This indicates that any 
silicate emission from the NLR must arise in its innermost regions. Alternatively, the silicate 
emission could come from the torus: we show that clumpy torus models give a reasonable 
fit to the silicate emission and the 2-20 /im SED, while at the same time obscuring the 
BLR. We emphasize that in the context of a clumpy torus the distinction of the transition 
from the outer torus to the inner NLR is more semantic than physical. Si multaneous mod- 



eling o f emission lines and silicate emission from the NLR, as proposed by ISchweitzer et al. 



(120081 ). may further illuminate the origin of the silicate emission features. Measurements on 



small spatial scales remain essential to identify the emitting structures in the cores of active 
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Fig. 3. — Clumpy torus models applied to the MIR spectra. Bars mark the regions employed 
in the fitting; IRS SL data at the short and long wavelength ends (black), and M i chelle 
measurements from 8-12.5 /im (gray). Photometric measurements ( lLawrence et al.1 11985 



Alonso-Herrero et al.lll996l . this work) and the IRS LL spectrum are plotted but not used 



in the fitting. The parameters of the best-fitting model 1 and model 2 are similar, with a 
thin torus (a = 15°), and limited optical depth through individual clouds (iy = 20), but 
the BLR remains obscured. Model 3 has a larger optical depth per cloud, a larger width 
parameter (cr = 45°), and fewer clouds along the torus equator. 




Fig. 4. — The clumpy torus can produce strong silicate emission even when the direct view 
of the AGN is blocked (model 4; i = 80°). This model yields a formally poor fit because the 
wavelength of the silicate features depends on the dust composition and grain size, which we 
do not vary. 
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galaxies. 
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